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Abstract 

A d-dimensional Ising model on a lattice torus is considered. As 
the size n of the lattice tends to infinity, a Poisson approximation 
is given for the distribution of the number of copies in the lattice of 
any given local configuration, provided the magnetic field a = a(n) 
tends to — oo and the pair potential b remains fixed. Using the Stein- 
Chen method, a bound is given for the total variation error in the 
ferromagnetic case. 
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1 Introduction 



The following situation, called "the law of small numbers" , is very classi- 
cal in probability theory. Suppose {/a}aeA is a finite family of indicator 
random variables, with the properties that the probabilities P{I\ — 1) are 
small and that there is not too much dependence between the Then, 
it is reasonable to expect the distribution of X)agA h to De approximately 
Poisson. In the theory of random graphs, inaugurated by Erdos and Renyi 
|S], such results are frequent (see 0] or PJ>i for a general reference). The 
I\ can, for instance, indicate the places in the random graph where a given 
subgraph appear. Some analogous results hold for random colorings of 
a lattice graph in dimension 2, corresponding to the context of random 
images ;6 . In both cases, the models are built on a large number of in- 
dependent random variables: the edges of a random graph or the pixels 
of a random image. In this article, we shall study Poisson approximations 
for sums of indicators defined from a large number of dependent random 
variables, namely the spins of an Ising model. 

Let us consider a lattice graph in dimension d > 1, with periodic boundary 
conditions (lattice torus). The vertex set is V n — {0, ...,n— l} d . The 
integer n will be called the size of the lattice. The edge set, denoted by 
E n , will be specified by defining the set of neighbors V(x) of a given vertex 
x: 

V(x) = {y^xeV n , \\y-x\\ p <p} , (1) 

where the substraction is taken componentwise modulo n, || • \\ p stands for 
the L p norm in R d (1 < p < oo), and p is a fixed parameter. For instance, 
the square lattice is obtained for p = p = 1. Replacing the L\ norm by the 
Loo norm adds the diagonals. From now on, all operations on vertices will 
be understood modulo n. In particular, each vertex of the lattice has the 
same number of neighbors. 

A configuration is a mapping from the vertex set V n to the state space W = 
{— 1, +1}. Their set is denoted by X n — W Vn and called the configuration 
set. Here, we shall deal with one of the simplest and most widely studied 
parametric families of random field distributions: the Ising model (see e.g. 

mns]). 

Definition 1.1 Let G n = (V n ,E n ) be an undirected graph structure with 
finite vertex set V n and edge set E n . Let a and b be two reals. The Ising 
model with parameters a and b is the probability measure p a ,b on X n = 
{— 1,+1} V ™ defined by: Ver G X n , 



Ha,b{o-) = ^— exp ( a ^2 a ( x ) + b ^2 °~{x)°~{y) ] , ( 2 ) 

ieV„ {x,y}£E n 



b 
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where the normalizing constant Z a ^ is such that Ma,&(°") = 1- 

Following the definition of |15| p. 2, the measure \i a .b defined above is 
a Gibbs measure associated to potentials a and b. Expectations relative to 
AVb will be denoted by IE a ^- 

In the classical presentation of statistical physics, the elements of X n are 
spin configurations; each vertex of V n is an atom whose spin is either pos- 
itive or negative. Here, we shall simply talk about positive or negative 
vertices instead of positive or negative spins. The parameters a and b are 
respectively the magnetic field and the pair potential. The model remain- 
ing unchanged by swapping positive and negative vertices and replacing a 
by —a, we chose to study only negative values of the magnetic field a. 
Various "laws of small numbers" have been already proved for the Ising 
model. Fernandez et al. EI] have established the asymptotic Poisson 
distribution of contours in the nearest-neighbor Ising model at low tem- 
perature (i.e. b large enough) and zero magnetic field. The Stein-Chen 
method is a useful way to get Poisson approximations; see pQ, [3J for a 
very complete reference or |Kj for the original paper of Chen. Barbour 
and Greenwood 2 have applied it to a class of Markov random fields; the 
bounds that they obtained for the Ising model are not quite explicit. In 
the same context as {HR, Ferrari and Picco ^I] have found bounds on the 
total variation distance between the law of large contours and a Poisson 
process. Ganesh et al. ^3] have studied the Ising model for positive values 
of b. Provided the magnetic field a tends to infinity, they proved that the 
distribution of the number of negative vertices is approximately Poisson. 
Our goal is to generalize the convergence in distribution given by ^Sj to 
any value of the pair potential b and to objects more elaborated than a 
single vertex. 

We are interested in the occurrences in the graph G n of a fixed local configu- 
ration r\ (see SectionElfor a precise definition and Figure^for an example). 
Such a configuration is called "local" in the sense that the vertex set on 
which it is defined is fixed and does not depend on n. Its number of occur- 
rences in G n is denoted by X n (rj). 

As the size n of the lattice tends to infinity, the potential a — a(n) will 
depend on n whereas the potential b will remain fixed. The case where 
a(n) tends to — oo corresponds to rare positive vertices among a majority 
of negative ones. As a consequence, the local configuration i] may occur 
or not in the graph, depending on its number of positive vertices k{r}). 
See Proposition 4.2 of for a precise description of this phenomenon. In 
particular, in order to get a nontrivial limiting result for the probability 
^aA x n(v) > 0), it is needed to take e 2a ^ of order n -d / fc M. Therefore, 
throughout this paper, the magnetic field a(n) will satisfy the identity 

e 2a(n) = cn- d l k ^ , (3) 
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where c is a positive constant. Our first result describes the asymptotic 
behavior of the number X n (rj) of occurrences of 77 in the lattice: it will be 
poissonian and will depend on k(rj) through @, but also on the geometry 
of 77 through its perimeter, denoted by 7(17). 

Theorem 1.2 Assume that the magnetic field a(n) satisfies (j3J and that 
the pair potential b is an arbitrary real number. As n tends to infinity, 
the distribution of X n (rj) converges weakly to the Poisson distribution with 
parameter c fc (»)) £-267(77) _ 

The proof is based on the moment method (see ^ p. 25 or Lemma 
I3.1|) . and requires estimates based on the local energy of 77 (Definition ^, lj l. 
The result of Ganesh et al. ^31 is obtained as a particular case when the 
pair potential b is positive and 77 is a single positive vertex: k(rj) — 1 and 
7(77) = 4. 

The Stein-Chen method makes it possible to obtain good estimates on the 
accuracy of Poisson approximations in terms of total variation distance. 
When the Gibbs measure fi a j, defined in satisfies the FKG inequality 
[T2*] (i.e. for positive values of the pair potential b), this method is applied 
to a sum of increasing random indicators (|24J) and produces Lemma |4.4I 
Then, bounds on the first two moments of the random variable X n (j]) 
(Lemmas 14.11 and 14. 511 allow to precise the Poisson approximation given 
by Theorem 11.21 This leads to our second result, where C{X) and V(X) 
respectively denote the distribution of X and the Poisson distribution with 
parameter A. 

Theorem 1.3 Assume that the magnetic field a(n) satisfies (0 and that 
the pair potential b is positive. Then, the total variation distance between 
C(X n {rf)) and the Poisson distribution with parameter e~ 2bl ^ satis- 
fies: 

The paper is organized as follows. The notion of local configuration 77 is 
defined in Section [21 Its number of positive vertices k{rj) and its perimeter 
7(77) are also introduced. Lemma l2"3l reduces proofs of Theorems ll.2l and ll.3l 
to clean local configurations. In this case, integers k{rj) and 7(77) naturally 
occur in the expression of the local energy of 77. Describing this quantity 
will be essential in our study. This allows us to control the conditional 
probability of 77 to occur in the graph (Lcmma l2.4|l . It immediatly follows 
that the expected number of occurrences of 77 in G n tends to c k ^ e~ 2hl ^ r, \ 
as n tends to infinity. Finally, Sections and 0] are respectively devoted to 
the proofs of Theorems 11.21 and 11.31 
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2 Conditional probability of a local configu- 
ration 



Let us start with some notations and definitions. Given a G X n = W n 
and V C V n , we denote by cry the natural projection of a over W v . If 
U and V are two disjoint subsets of V n then o~u<?'v is the configuration on 
U U V which is equal to a on U and a' on V. Let us denote by 8V the 
neighborhood of V (corresponding to Q): 

SV = {ye V n \ V, 3x e V, {x,y} e E n } , 

and by V the union of the two disjoint sets V and 8V. Moreover, \V\ 
denotes the cardinality of V and T{V) the cr-algebra generated by the 
configurations of W v . 

As usual, the graph distance dist is defined as the minimal length of a path 
between two vertices. We shall denote by B(x,r) the ball of center x and 
radius r: 

B(x, r) = { y e V n ; dist(x, y) <r} . 

In the case of balls, B(x,r) = B(x,r + 1). In order to avoid unpleasant 
situations, like self-overlapping balls, we will always assume that n > 2pr. 
If n and n' are both larger than 2/?r, the balls B(x, r) in G n and G n > are 
isomorphic. Two properties of the balls B(x,r) will be crucial in what 
follows. The first one is that two balls with the same radius are translates 
of each other: 

B(x + y, ?') = y + B(x, r) . 

The second one is that for n > 2pr, the cardinality of B{x 1 r) depends only 
on r and neither on x nor on n: it will be denoted by f3(r). The same is 
true for the number of edges {y, z} S E n with y, z 6 B(x, r), which will be 
denoted by a{r). 

Let r be a positive real, and consider a fixed ball with radius r, say 
B(0,r). We denote by V r = W B ^ ' rS> the set of configurations on that 
ball. Elements of T> r will be called local configurations of radius r. A local 
configuration rj e V r is determined by its subset V+(ry) C B(0, r) of positive 
vertices: 

V-+(77) = {xeS(0,r), V (x) = +l} . 

The cardinality of this set will be denoted by k(r/) and its complementary 
set in B(0,r), i.e. the set of negative vertices of ry, by V-(r]). Of course, 
there exists only a finite number of local configurations of radius r (precisely 
2/ 3 ( r )). In what follows, 77, 7/ will denote local configurations of radius r 
and C' those of radius larger than r. 

A local configuration 77 6 D r is said clean if its subset of plus vertices V+ (77) 
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is included in the ball B(0, r — 1). In other words, vertices of a clean local 
configuration which are at distance r from the center are negative. Figure 
n shows an example of such a local configuration. 



----- + - 
--:+-:++- 
--- + + -- 
- + - - - + - 
_+__:+__ 



Figure 1: A clean local configuration r\ with k(rf) = |V+(r7)| — 10 positive 
vertices, in dimension d — 2 and on a ball of radius r = 3 (with p = 1 and 
relative to |j • ||oc). 

Let rj G T> r . For each vertex x G V n , denote by r\ x the translation of rj 
onto the ball B(x,r) (up to periodic boundary conditions): 

Vy G V n , dist(0, y) < r r\ x (x + y) = rj(y) . 

Let us denote by IJj. the indicator function defined on X n as follows: ZJ?(cr) 
is 1 if the restriction of the configuration a G X n to the ball B(x,r) is r) x 
and otherwise. Finally, let us define the random variable X n (r)) which 
counts the number of copies of the local configuration rj in G n : 

Due to periodicity, this sum consists of n d indicator functions 1%, which 
have the same distribution. 

In order to control the random variable X n (rj), we describe its "local be- 
havior" by introducing the local energy of 77. Let us start with the following 
definition. 

Definition 2.1 Let x G V n and a G W B(x ' r+1 'K The local energy H B{ - X ^(a) 
of the configuration a on the ball B(x,r) is defined by: 

H B ^ r \a) = a(n) £ a(y) + b £ a(y)a(z) , 

y£B(x,r) {y,z}£E n 

(y£B(x,r))V(z£B(x,r)) 
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where (y £ B(x,r)) V (z £ B(x,r)) means at least one of the two ver- 
tices y and z belongs to B{x, r) (the other might belong to its neighborhood 
5B{x,r)). 

Let us fix a vertex x and denote merely by B the ball B(x, r). For any 
local configuration 77 £ V r and for any a £ W SB , the local energy H B (r] x a) 
on B of the configuration which is equal to rj x on B and a on SB can be 
expressed as: 



a(n)(2k(r])-l3(r)) + b 



( _ \ 

Vx(y)Vx(z) + ij x (y)a(z) 

\ y,z£B y£B,z£SB / 



E 



(4) 



Actually, this notion of local energy allows us to explicitly write the con- 
ditional probability H a ,b{^S = a e W SB : 



Ha,b(Z> = Ik) = ^ pg B (<er) ■ (5) 



As we shall see in Lemma 12.41 bounding the above conditional probability 
is central in our study. 

An easy way to connect the number of copies of the local configuration 77 
to its local energy consists in writing, for any given vertex x: 

Ea, b [X n (ji)] = E a , b [n d I2] 

= E a , b [n d ^ a . b {P x = l\T{5B))] . (6) 

Here, \x a , b {l x — M-F(SB)) represents a JF(<5B)-measurable random vari- 
able and, for a £ W 5B , fi a , b (l% = 1| T{5B)){&) = H a ,b(IS = l W) a con- 
ditional probability. Note that the set SB has bounded cardinality (not 
depending on n) . Then, from a convergence result on the random variable 
nd f-a,b(IS — 1\3~(5B)) it wm be easy to obtain a similar result for its ex- 
pectation, i.e. for lE a , b [X n (rj)]. 

We will now give the reason for the hypothesis (JHJ) , that links the mag- 
netic field a(n) to the number of positive vertices of the local configuration 
77. The event X n (77) > corresponds to the appearance of 77 in the graph 
G n . In [7j Proposition 4.2, it has been proved that: 

if lim e Mn)k( v ) n d _ then 

71 — >00 

lim Ha, b (Xn(v) > 0) - and lim JE a . b [X n ( V )} = ; (7) 
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if lim e Mn)H v ) n d = +00 then 

n-^-oo 

lim (i a , b (X n (rj) > 0) = 1 and lim E a , b [X n (r})] = +00 . (8) 

n — >oo n — >oc 

In particular, the element of T> r having only negative vertices, called the 
null local configuration and denoted by rf , has a probability which always 
tends to 1. 

From now on, assume that 77 has at least one positive vertex; k(rj) > 1. 
Using the vocabulary of the random graph theory, statements (Q) and (JHJ 
essentially mean that the quantity 

n -d/k(v) ig 

the threshold function (for 
e 2a(n)^ Q £ p r0 p er ty X n {rf) > 0. It does not depend on the radius r of 
the ball on which the local configuration 77 is defined: r is just a phan- 
tom parameter which serves only to ensure that 77 is a local configuration. 
Actually, the function n~ d ' k ^ Tl ' only depends on the number of positive 
vertices of 77. Roughly speaking, if e 2a (") is small compared to n~ d ^ k ^ v \ 
then asymptotically, there is no copy of r\ in G n . If e 2a (") is large compared 
to n~ d l k ( r, \ then at least one copy of 77 can be found in the graph, with 
probability tending to 1. 

Consequently, in order to get a nontrivial limiting result for the probability 
(J>a,b(X n (r)) > 0), it is needed to take e 2a< ^™^ of order n ~ d / k ^'> . Hence, for 
the rest of this article, is satisfied, i.e. 

e 2a(n) _ cn -d/k( n ) 

for some positive constant c. Under this hypothesis, statement JJJ says that 
asymptotically there will be no local configurations with (strictly) more 
than fc(ry) positive vertices in the lattice. The following lemma quantifies 
this result. 

Lemma 2.2 Let 77 6 T> r and suppose that e 2a (") = cn~ d / k ( v \ for some 
constant c > 0. Let R > r an integer and £ G "Dr. Then, there exists 
a constant M\ > such that for all n, for all vertex x € V n and for all 
configuration a G X n , 

nV,b(4 = lkiB(x,fl)) < M ie 2a (")W«- fc W) , (9) 

and for all n, 

E a ,b[ X n(0] < Mie 2a ^ k ^~ k ^ . (10) 

Proof. Let x G V n and denote merely by Br the ball B(x,R). For any 
configuration a G W SBr , the local energy H BR (( x a) is given by 10}. Since 
the set Br has a bounded cardinality, there exists a constant K > 0, only 
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depending on the pair potential b and the radius R (and not on n, nor x, 
nor a), such that 

H B *iQ x a) - H B *(Ca) < 2a(n)k(() + K , 

where C° is the null local configuration of radius R (k(C,°) = 0). Hence, by 
relation J5J and for all a G W SBlt : 

< e 2a(n)fc(C)+K _ 

Finally, hypothesis (0 provides the first inequality of Lemma \l. 21 

with Mi = c k ^e K > 0. The quantity n d fj, a ^(I^, = l\a) is bounded uni- 
formly on the configuration a e W SBr . So, its expectation satisfies the 
same inequality and 11U|) follows. □ 



A primary consequence of Lemma 12.21 consists in reducing our study 
to clean local configurations. To any given 77 £ XV, a local configuration 
77 £ £V+i is associated by the following process: 



f]{x) 



T](x) if x e 5(0, r) , 
-1 if dist(x,0) = r + 1 



The local configuration 77 is clean and satisfies fc(7/) = fc(?y). Note that the 
inequality X n (rj) < X n {r)) holds for all size n. Actually, these two random 
variables are asymptotically equal. Indeed, assume that r\ occurs on the 
ball B(x, r). Then, hypothesis J3J forces vertices at distance r + 1 from x 
to be negative with probability tending to 1: 

lim Ma , 6 (I* = l\F x = 1) = 1 . 

Lemma 12.31 expresses this result in terms of total variation distance. Recall 
that if /_t and v are two probability distributions, the total variation distance 
between \i and v is 

d T y(/-«, v) = sup \fi(A) - v{A)\ , 

A 

where the supremum is taken over all measurable sets. Besides, the prob- 
ability distribution of a random variable X will be denoted by C{X). 
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Lemma 2.3 Let rj G T> r and suppose that e 2a ^ = anr d ' k ^ n ' , for some 
constant c > 0. Then, the total variation distance between distributions of 
X n (fj) and X n (j]) satisfies: 

d TV [C{X n {rj)),C{X n {ri))) = 0(n~ d / k ^) . (11) 

Furthermore, the difference between the expected numbers of copies of local 
configurations rj and rj in the graph G n tends to 0: 

lim \lE a . b [X n (ij)} - E afi [X n {n)]\ = . (12) 

n — >+oo 

In this paper, if f(n) and g(n) are two positive functions, notation 
f(n) = 0(g(n)) means that there exists a constant C > such that, for all 
n, f(n) < Cg(n). 

Proof: On the one hand, let us introduce the subset T)]l +1 of V r+ i defined 
by: 

= K G V r+1 , Vz G B(0,r), ((x) = V (x)} . 

The local configuration 77 is the only element of £>J? +1 satisfying k(fj) = k(j]), 
all the others having at least k(rf) + 1 positive vertices. Moreover, the sum 
of all copies of elements of V^, +1 is equal to the number of copies of 77: 

x n {r,) = x n (f,) + J2 ■ ( 13 ) 

C6X»? +1 \{17} 

On the other hand, the total variation distance between two probability 
distributions can be written as 

d-TV (Mi v) = inf{F(X ^ Y), C(X) = fj, and C(Y) = v} . (14) 

Using this characterization and the identity $Tj$, it follows: 

(hv{C{X n {fi)),C{X n {rj))) < fi a , b (X n (fi) ? X n (r])) 

< fi a ,b(X n (v) > X n {fj)) 

< ^ b J2 *»(0>o] . 

\Cex>;' + i\{>j} / 

The above sum is an integer valued variable. So, its probability of being 
positive is bounded by its expectation. Hence, 

d TV (C(X n (fi)),C{X n (r)))) < ®a,b[X n {0] 

Cex>;> +1 \{»)} 

< \V ,1 r+l \M 1 e 2a ^ , 
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by Lemma l2~2l Thus, using e 2a ^ = cn d / fe (*w, we deduce that the to- 
tal variation distance between the distributions of X n (fj) and X n (j]) is a 
0( n -«V*(i)). 

Finally, the same is true for the absolute value of the difference between 
the expectations of X n (fj) and X n {rf) since: 

\lE atb [X n (fi)} - M a<b [X n {ri)]\ = JEa,Mn{Q] • 

Ce^ +1 \«} 

□ 

If the random variable X n {f]) converges weakly as n tends to infinity to 
a limit v 1 inequality implies that the same is true for X n (rj). Conse- 
quently, replacing r with r+1 and without loss of generality, we can assume 
that vertices of the reference ball B(Q,r) which are at distance r from the 
center 0, all belong to V-(rf) (as in Figure^]), i.e. r\ can be assumed clean. 

Now, the geometry (in the sense of the graph structure) of the set V+(rj) 
of positive vertices of 77 takes place in our study. Precisely, let us define the 
perimeter 7(77) of a local configuration rj G T> r by the formula: 

l(v)=V\V+( V )\-2\{{x,y}eV + (r 1 )xV + (r l ), {x,y} G E n }\ , 

where V represents the number of neighbors of a vertex. In particular, the 
perimeter of a local configuration is always an even integer. For instance, 
that of Figure Q is equal to 58. If 77 is clean, its perimeter 7(77) merely 
becomes: 

7(»0 = \{{x,y} e ^+(77) x V1(t7) , {x,y} G E n }\ . 

In this case, 7(77) represents the number of pairs of neighboring vertices x 
and y of B(0, r) having opposite spins under 77. As we shall see in the proof 
of Lemma 12.41 the perimeter of a clean local configuration easily occurs in 
the expression of its local energy. It is the reason why we reduce our study 
to that of clean local configurations. 

The following lemma will play an essential role in the proofs of Theorems 
11.21 and ll. 31 it gives a uniform bound for the random variable n d /i a .f,(/J = 
l\F{6B(x,r))). 

Lemma 2.4 Let rj be a clean local configuration of radius r and suppose 
that e 2a ^ — cn~ d / k ^ , for some constant c > 0. Then, there exists a 
constant M2 > such that for all n, for all vertex x G V n and for all 
configuration a G X n , 

c fc We -2 67(l) ) (1 _ M2e 2a(»)j < n d^ b (jr, = l\ GsB{jB r) ) < c Hv) e -2bj( V ) > 

(15) 
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and for all n, 

C fc(r,) e -26 7 ( I) )( 1 _ M 2 e 2a ^) < JE atb [X n (rj)} < c^e" 26 ^ . (16) 

Since the quantity e 2a (") tends to as n tends to infinity, the inequalities 
(|15fl and HKifl yield the two following limits. For any vertex x and any 
configuration a, 

lim n d fx a , b (I2 = l\<r 5B(x , r) ) = lim JE a , b [X n ( V )] = c k M e ~ 2b ^ . (17) 

Thanks to l|12(l of Lemma 12.31 the latter limit is valid for any element of 
T> r (not necessary clean). 

Proof: Let a; be a vertex of V n and denote by B the ball B(x, r). Since the 
expectations of the variables X n (n) and n d \x a ^ij x = 1\T{6B)) are equal 
(see ©), (|16|l is an immediate consequence of Ijl5|) . So, let us prove this 
relation. 

Let us start with inserting the perimeter 7(77) in the expression of the local 
energy of 7/. Assume that r\ x occurs on B. Then, there are 7(77) edges 
{y,z} e E n with y,z e B satisfying r] x (y)ri x (z) = -1 and a(r) - 7(77) 
ones satisfying rj x (y)i] x (z) — 1. Hence, for all a e X n , the local energy 
H B (t] x o-5b) can be expressed as: 



a(n)(2k(r))- (3(r)) + b 



( \ 

a(r)- 27(77)+ (-!>«»(*) 

\ y£B,zeSB / 



The factor (—1) in the latter sum comes from the fact that, by hypothesis, 
vertices at distance r from x are all negative. Let 77' G T> r be a local 
configuration of radius r with A; (7/) positive vertices. Then, the difference 
H B (v' x °~sb) — H B (r] x asB) between the local energies of 77^. and 77^ is equal 
to: 

H b ( v ' x <j sb ) - H B ( Vx a SB ) = 2a(n)(fc(V) - Hv)) + b (27(77) + Q( V ' X )) , 
where 

Q(V'J = vL(yW x ( z ) ~a(r) + Y {rfxiv) + t)<T5B(z) . 

{y,z}eE n { s ,4eB~ 

y,z£B y£B,z£5B 

The real 2(77^,) does not depend on n and can be bounded uniformly on 
configurations 77^., a: |S(t7^.)| < 2a(r+l). Moreover, note that the null local 
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configuration of radius r satisfies Q(r] x ) = 0- So, using (j^J, the quantity 
H B {if x asB) - H B {rj x a SB ) becomes: 

H B {rP x a SB )-H B ( WB ) = -2a(n)k(r 1 ) + 2b 1 (j 1 ) 

= lo s(^))+ 2 ^), 

for any configuration a £ X n . Then, using the explicit formula for the 
conditional probability /u Q , — ^\&8b) (relation 10)), we get: 

HM = ^8B) < e» B ^B )-» B (vl« SB ) 

c fe(r,) e -26 7 (> ) ) 
< ; . 



i.e. the upper bound of l|15|) . The lower bound is obtained as follows. For 
any configuration a 6 X n : 



n d 



V P H B (v>SB)-HB( V;c a SB ) 
c fe(77) e -2bj(r,) 

± T Z-iri' ,k(ri>)>0 c 

> 



c k(ri) e -2b~,(ri) 



l + e 2o(n) | X) r | e 2|6|a(r+l) ' 

Let Mi denote the quantity |r> r |e 2 ' b ' Q<r+1 " ) ; it only depends on the pair 
potential b and the radius r. Finally, the inequality 

Vu > -1, — !— > 1 - u 
1 + u ~ 

implies the lower bound of (|T5|) . □ 



3 Poisson approximation 

This section is devoted to the proof of Theorem 11.21 the distribution 
of X n {rf) converges weakly to the Poisson distribution with parameter 
c M'i)e- 21 Tw), Previous notations and hypotheses still hold. In particu- 
lar, the relation J3J) between the magnetic field a(n) and the number k(rf) 
of positive vertices in r\. 

Before proving this result, it is worth pointing out here the role of the pair 
potential b. First, remark that local configurations of radius r having the 



13 



same number of positive vertices can have different perimeters. Theorem 
11.21 assures that the probability (for fi a ,b) of the local configuration 77 of 
occurring in the graph is asymptotically equal to 



1 -e 



_ c kM e -2b-r(v) 



So, if b > (resp. b < 0), this asymptotic probability is a decreasing (resp. 
increasing) function of the perimeter 7(77). In other words, if b > (resp. 
b < 0), among the local configurations having the same number of positive 
vertices, those having the highest asymptotic probability of occurring in 
the infinite graph are those having the smallest (resp. largest) perimeter. 
If the pair potential b is null then the perimeter 7(77) of the local configu- 
ration 77 has no influence. All local configurations having the same number 
of positive vertices have the same asymptotic probability 1 — e~ c of oc- 
curring in the graph. In the 2-dimensional case, this is Theorem 2.4 of 

In order to prove Theorem 11.21 we use the moment method based on 
the following lemma (4 p. 25). 

Lemma 3.1 Let (Y n ) ne fq* be a sequence of integer valued, nonnegative 
random variables and X be a strictly positive real. For all n, I £ N* define 
Mi(Y n ), the l-th moment ofY n , by 



If, for all I £ N* , liniyj^oo Mi(Y n ) = X 1 then the distribution ofY n converges 
weakly as n tends to infinity to the Poisson distribution with parameter X. 

First, note that Lemma l2~3l reduces the proof of Theorem ll.2l to a clean 
local configuration rj. 

So as to lighten formulas, the quantity c k ^ e~ 2bl ^ will be simply de- 
noted by A. Thanks to Lemma 13.11 we just need to prove the conver- 
gence of Mi(X n (rf)) to X 1 for every 1 £ N*. The case I = 1 has already 
been treated at the end of the previous section, where it was proved that 
Mi(X n (r])) = lE a fi[X n (ri)] tends to A. From now on, fix an integer I > 2. 
In our case, the variable X n (n) counts the number of copies in the graph 
G n of the local configuration rj. Then, the quantity Mi(X n (rf)) can be in- 
terpreted as the expected number of ordered Z-tuples of copies of r/. 
If B represents a set of balls of radius r whose centers belong to V n , then 
two elements B(y, r) and B{z, r) of B will be said to be connected if there 
exists an integer m and balls B%, . . . , B m £ B such that B\ = B{y,r), 



Mi{Y n ) 



E[Y n {Y r 



n 



l)...(Y n -l + 1)] 




14 



B m = B(z, r) and for j = 1, . . . , m— 1, Bj n Bj+i ^ 0. This last condition 
allows balls Bj and Bj+i to be disjoint but their centers are at distance 
from each other at most 2r + 1 . The connectivity is an equivalence relation 
on the set B. 

For s = 1, . . . , I, denote by Cz(s) the set of Z-tuples of vertices (sci, . . . ,x{) 
such that the set of balls {B(x\,r), . . . , B(xi,r)} is composed of s equiva- 
lence classes for the connectivity relation. Then, the l-th moment of X n (r]) 
becomes: 



(xi,...,Kj)ec ( (s) 



II x ■ ■ ■ x II 



The term corresponding to s — I in the above sum will be denoted by 
M[(X n (r])) and the remaining sum by M ; (X n (rf)). We are going to prove 
the two following limits 



lim M[(X n (r])) = \ l 



(18) 



lim M l (X n (r))) = , 

n — *oo 

from which Theorem 11.21 follows . 

Let us first check the cardinality of Ci(s). 

Lemma 3.2 There exists a constant C > such that 

Vs = l,...,l, \Ci(s)\ <Cn ds . 
Furthermore, the cardinality of Ci(l) is equivalent to n dl : 

\Ci(l)\ 



lim 



ill 



1 



(19) 



(20) 



(21) 



Proof. Let (x\, . . . , x rn ), m < I, be a m-tuple of vertices such that the set of 
balls {B(xi, r), . . . , B{x m , r)} is composed of only one equivalence class for 
the connectivity relation. Each vertex Xj, j = 1, . . . , m, necessary belongs 
to the ball B(xi,R) with R = l(2r + 1). So, the number of such m-tuples 
(afi, ... ,x m ) is bounded above by (3{R) l ~ 1 n d . Applying this argument to 
each of the s equivalence classes of the set formed of the I balls centered at 
vertices of a given /-tuple {x\, . . . ,xi) £ C/(s) provides: 

\Ci{s)\<f3{R) s ^ 1) n ds . 

Note that C = 0(R) 1 ^ 1 -^ only depends on integers / and r and satisfies 
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Now, we want to choose I vertices x%, . . . , xi such that (x\, . . . , xi) € Ci(l). 
For the first vertex xi, there are n d possibilities. Let 2 < j < I and 
suppose vertices Xi, . . . ,Xj-i have been chosen. For the j — th choice, the 
set of all vertices x such that dist(x,Xk) < 2(r + 1) for some 1 < k < 
j — 1, must be avoided. The cardinality of this set is bounded above by 
(j — 1) x /3(2(r + 1)) whatever xi, . . . , Xj-i. This bound does not depend 
on n. As a consequence, from the inequalities 

[] (n d -(3-l)P(2r + l)) < \C,(l)\ < n dl , 

i<j<i 

estimate l|2"T|l follows. □ 



Let us prove the hrst limit (|18|l . Let (xi, . . . ,xi) £ Ci(l). By definition 
of the connectivity relation, note that, for 1 < i, j < I and i j, no vertex 
of the ball B(xi,r) can be a neighbor of a vertex of the ball B(xj,r). 
The Gibbs measure fi a ,b yields a Markov random field with respect to 
neighborhoods defined in Q (see for example 15 , Lemma 3 p. 7). As a 
consequence, 



Ma, 6 



\i=i 



En. 



IE,, 



En. 



1 



\i=l 



F(u l j=1 8B( Xj ,r)) 



H^ a , b (ll = l\^(5B( Xi ,r))) 



(22) 



Since the local configuration rj is clean, Lemma l2~^l and relation ill'L'l) provide 
a control of the probability ^ a ,b{Ix 1 x ■ • ■ x 13, — 1) : 



1 < 



uniformly on the /-tuple {x\, . . . , x{) 6 Ci(l). Hence, 



\Ci(l)\ 

n dl 



X 1 (l - M 2 e^) 1 < M((X n (r,)) < l -^X l . 



Finally, as the ratio \Ci(l)\ divided by n dl tends to 1 (relation H21J) L the 
quantity M[(X n (ri)) tends to the searched limit. 



1G 



There remains to prove that M l (X n (rf)) tends to as n tends to infinity. 
The intuition is that if the local configuration rj occurs on two balls B(x, r) 
and B(x',r) with dist(x,x') < 2(r + 1), then locally (strictly) more than 
k{rj) positive vertices are present in a ball of radius 2(r + 1). This has 
vanishing probability by Lemma 12.21 

Let us prove that every term of the sum defining M ; (X n (rf)) tends to 0: 
fix an integer 1 < s < I — 1. Let {x\, . . . , x{) £ C;(s). The set of balls with 
radius r, centered at these vertices, splits into s equivalence classes, say 
EC(1), . . . , EC(s). Let us denote by Cj the union of balls belonging to the 
equivalence class EC(j). Once again, we use the markovian character of 
the Gibbs measure [i a y- 



Va,b 



n^=i 



T[»a,b\ n T l = 1 

3 = 1 \i,B(x,,r)£EC(j) 



As a consequence of s < l—l, there exists at least one connected component, 
say -EC(l), having at least two elements: let x{l) and x'(l) be two vertices 
satisfying dist(x(l), x'(l)) < 2r + 1 and B{x{l),r), B{x'{l),r) S EC{1). 
For every j = 2, . . . , s, denote by x(j) one of centers of balls belonging to 
EC(j). Then, we can write: 



i,B(x it r)eEC(j) 



< JEa,b 
X 

~ n d 



U) 

IE a ,b 
JEa.b 



T{SC 3 U 6B(x(j), r)) T{5C 3 ) 
F(5B(x(j),r))} F(5Cj) 
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by Lemma 12.41 (n is clean). This last inequality allows us to write: 



eEC(i) 



n 



. i,B(xi,r) 
\ GBC(l) 



n ^=1 



•F(tfCi) 



J 



< 



d fJ-aA 1 ^!) - ^'(1) _ 1 ) 



e) set of local configurations 
having at least k(-q) + 1 positive vertices. Then, the event 1%,^ = I^m = 1 



(23) 



Denote by V^+i the (finite) set of local configurations of radius 2r + 1 



implies that one of the elements of 'D^r+i occurs in B(x(l),2r + 1). It 
follows that: 



c v >kM 



< n- d \V> r k _l?\ M ie 2a ^ , 

by Lemma l2~2l As a consequence, the following bound does not depend on 
the ^-tuple (x\, . .. ,xi) GCi(s): 



Finally, Lemma 13.21 implies that 

lEa.b 



X...XII 



KCX*- 1 |X>OT| Ma e 2Q (") 



(ai,...,x;)eC( (s) 

which tends to as n tends to infinity. Theorem 1 1 . 21 follows . 



4 The ferromagnetic case 

In this section, we suppose the pair potential b is nonnegative and the 
magnetic field a(n) satisfies the relation Q: 

e 2a(n) _ cn -d/k{n) 
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for some positive constant c > 0. Under these hypotheses, Theorem 11.31 
says that the total variation distance between C{X n {rj)) and its Poisson 
approximation V{c k ^e-^^) is a 0(n~ d / k ^). 

We believe that n~ d ^ k ^ is the real speed at which the total variation 
distance between C{X n (r])) and V(c k ^e- 2b ^) tends to zero. Indeed, it 
seems to be true for the upper bound given by Lemma Em ffor more details, 
see Chapter 3 of 3 ). Besides, in the case where the local configuration 77 
represents a single positive vertex (with £;(?/) = 1, 7(77) = 4, p = 1 and 
relative to || • ||i), Ganesh et al. ^3] proved that 

logn~ d 

as n tend to infinity. 

The rest of this section is devoted to the proof of Theorem ll.3l First, 
note that Lemma \'2 .31 reduces the proof to a clean local configuration 77. 
Let us start with some notations and definitions. There is a natural partial 
ordering on the configuration set X n = { — l,+l} v ™ defined by a < a' if 
o(x) < c'(x) for all vertices x € V n . A function / : X n — > 1R is increasing 
if /(c) < /0 7 ') whenever a < a' . 

From the local configuration 77, let us define the subset T> r (r]) of V r by: 

25,(77) = {7/ e2? r ,U+(r/) 3^(77)} . 

Each local configuration of 2?*(?7) = £V(?7)\{?/} has at least ^(77) + 1 positive 
vertices. Moreover, by definition of 25,(77) and for all x G V n , the indicator 
defined by 

T"= V 

X / j X 

ri'eVrin) 

is an increasing function. Let us introduce the corresponding random vari- 
able X n (r]): 

Xn(ri) = E 7 " ( 24 ) 

xev n 

= x n ( v ) + J2 x »tf) > ( 25 ) 

whose expectation IE a ^[X n (j])] will be simply denoted by A„. As in the 
previous Section, the quantity c fe (' 7 - ) e -2b7 '-''' 1 will be denoted by A. 
The proof of Theorem 11.31 is organized as follows. The total variation 
distance between C(X n (rj)) and V(X) is bounded by: 

d T v{C{X n (r,)), C{X n {ri))) + d TV {C(X n {ri)),T{\ n )) + d TV {V{\ n ), V{\)) . 
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Let us respectively denote by T\, T2 and T3 the three terms of the above 
sum. We are going to prove that each of them is of order 0{n~ d / k ^). 
Terms T\ and T3 are respectively dealt with using Lemmas 14.11 and 14.21 
Applied to the family of indicators {fl,x 6 V n }, the Stein-Chen method 
gives an upper bound for T2 fLemm.a l4.4jl . Finally, Lemma l4~5l implies that 
this upper bound is a 0(n~ d ^ k ^). 

Hypothesis © implies that occurrences of local configurations with 
(strictly) more than k{rf) positive vertices have vanishing probability. Hence, 
the random variables X n (r]) and X n (rj) will be asymptotically equal. So 
do their expectations. 

Lemma 4.1 The total variation distance between the distributions of X n (rj) 
and X n {ri) satisfies: 

d TV (C(X n (r,)),C(X n (r,))) = (D(n~ d / k ^) . (26) 

Furthermore, there exists a constant M3 > such that for all n: 

A(l - A/ 3 e 2a (")) < A„ < A(l + M 3 e 2a{n) ) . (27) 



Proof. Thanks to relation between X n {rf) and X n {rj) l|25|l and characteri- 
zation l|14(l of the total variation distance, we get: 

d TV (C(X n {r))), C(X n (r)))) < fi a ,b(Xn(v) ^ X n {r,)) 

< fi a>b X ^') > 

< \V*(rj)\ M ie 2a ^ , 

by Lemma l2~2l So, the quantity dxv {£(X n (v))> £(X n {v))) is equal to 

Using the previous inequalities, a control of the expectation A„ of X n (rj) 
is obtained: 

1E a , b [X n {n)] < A„ < E atb [X n (r,)} + |P r *(77)|M ie 2a (™) . 
Since n is clean, relation (|16|l can be applied. The above control becomes: 

A (l - M 2 e 2Q(,l) ) < A„ < A (l + \V* r (r 1 )\M 1 \- 1 e 2a ^ . 
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Recall that constants Mi and Mi do not depend on the size n. Relation 
(|2*7jl follows by letting: 

M 3 = max{Af 2 , ^(^MiA" 1 } . 

□ 



Using (|27|l . we shall now bound T3. 

Lemma 4.2 TTie £o£aZ variation distance between the Poisson distribu- 
tions with parameters X n and X satisfies: 

T 3 = d TV (V(X n ),V(X)) = 0{n- d ' k ^) . 

Proof. The total variation distance between two probability distributions 
on the set of integers can be expressed as: 



d TV (fi, v) = - ImM - v{i 



2 

m>l 

Let to > 1. Thanks to relation f2T| . the difference A™e _A ™ — X m e~ x is easily 
controlled. Thus, dxvCPi^n), T'(A)) is bounded by |E m >i max { a 'n>^'n} 
where: 

am = f e AA/ 3 e-<") (1 + M3e 2a(n) )m _ ^ 

to! V / 

and 

(} m = (l - e -AM 3 e-'") (1 _ M3e 2a(n))^ . 

m! V / 
Using the convexity of the function 

f m ■ } - 1, 1[— > M, x .— » (1 + z) ro e A * , 

one easily checks that a m > /3 m , for all m £ IN*. As a consequence, 

m>l 



which is of order 0(n d / k M) by relation 10. □ 

There remains to bound the term T 2 = dTv{£{X n {rj)), V{X n )). This is 
based on the Stein-Chen method and particulary on Corollary 2.C.4, p. 26 
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of P] which is described below (Proposition 14. 3|1 . Let be a family 

of random indicators with expectations Tti . Let us denote 

W = J2 T i and 9 = J2*i- 

i£l iel 

The random variables are positively related if for each i, there exists 

random variables {Jj,i}jei defined on the same probability space such that 

£(Jj,i, j 6/)=% je/|I< = l) 

and, for all j ^ i, Jj.i > i}. 

Proposition 4.3 // £/ie random variables are positively related 

then: 

1 _ P -e ( 

d T v(£{W),V{6)) < \Var(W)-9 + 2j2^ 

Proposition 14.31 can be applied to our context. Indeed, for a positive 
value of the pair potential b, the Gibbs measure ^, a ,b defined by 10) satisfies 
the FKG inequality, i.e. 

lE a ,b[f9] > lE a , b [f}E a , b [g} , (28) 

for all increasing functions / and g on X n : see for instance Section 3 of 
[EJ. Then, Theorem 2.G p. 29 of implies that the increasing random 
indicators {T^, x £ V„} are positively related. Replacing Ii with J^, W 
with I„(ij) and 9 with A„, Proposition 14.31 produces the following result. 
This is the only place where the hypothesis b > is actually used in the 
proof. 

Lemma 4.4 // the pair potential b is nonnegative then the following in- 
equality holds: 

d T v(C(X n (r,)),V(X n )) < ( Var a 4X n ( V )] - A„ + 2 V dE a ^ x f ) . 

Xn \ x^L J 

(29) 

The bound (|29|l indicates that, as in +oo, the sum J2 x gv ^a.b[I x ] 2 
is small and the distance to the Poisson approximation is essentially the 
difference between the variance and the expectation of X n (ij) . Using good 
estimates on the first two moments of the random variable X n (r]), this 
difference will be bounded. The case of the first moment of X n (i]), i.e. 
its expectation A„, has been treated in Lemma |4.1I The following result 
concerns its second moment: 
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Lemma 4.5 The second moment M 2 (X n (i~i)) — lE a ^[X n (rf)(X n {rj) — 1)] 
of the random variable X n (ij) satisfies: 

M 2 {X n {i 1 )) = \ 2 + 0{n- d ' k ^) . 

Writing the variance of the variable X n (rj) as M 2 (X n (rj)) + \ n — A 2 and 
the sum J2 x gv ^a,,b[fl] 2 as the ratio X 2 l /n d , we deduce from Lemma IPI 
that: 

1 



T 2 < — lM 2 (X n ( V )) + \i(— -1) 

The inequalities given by Ij27(l and Lemma f4.5l allow us to control the ex- 
pectation A„ and the second moment M 2 (X n (i])) of the random variable 
X n (i]). This implies: 

T 2 = 0(e 2a ^) + ^ , 
n a 

which is a 0(e 2a ^) = 0{n- d ' k ^) since k(r]) > 1. 

Let us finish the proof of Theorem II. 31 by proving Lemma 14.51 

Proof (of Lemma \4-5\) : First, recall that C 2 (s), for s = 1,2, represents the 
set of couples (xi,x 2 ) whose set {B(xi, r), B(x 2 , r)} splits into s equiva- 
lence classes for the connectivity relation. In other words, (xi,x 2 ) belongs 
to C 2 (l) if dist(xi,x 2 ) < 2r + 1 and to C 2 (2) otherwise. So, the second 
moment of X n (r)) is equal to: 



M 2 (X n {r{)) = Y,JEa,b 



E 

(x 1 ,x 2 )&C 2 (s) 



1 XI X 



Each indicator 7^ is defined as the sum of FJ. , rj' £ T> r (rj). Hence, the 
second moment M 2 (X n (rj)) becomes: 



M 2 (X n ( V )) = (E 1 (r n ,r ]2 )+E 2 (7 ll , m )) 
where for s = 1,2, the quantity E s {rji^rj 2 ) is defined by: 



(30) 



E s { m ,7 l2 ) = m a 



E 



r'fi jm 



(xi,X 2 )&C 2 (s) 



Let (771,772) be a couple of local configurations belonging to T) r (rj) and 
(xi,X2) be a couple of vertices. In a first time, consider (xi,x 2 ) € C 2 {1). 
Then it has been already seen at the end of the previous section that the 
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event = = 1 implies that one of the elements of ^2r+i necessary 
occurs in B(xi, 2r + 1). It follows that: 



< n- d Y, m »AXn{Q] 
- d \V> r k lf\M ie 2a ^ , 



< n 



by Lemma l2~2*l Thus, thanks to Lemma l3~2l we deduce that £^(771,772), for 
s = 1, 2, is a 0(e 2a (™)) = 0(n- d / fc ^). 

Now, let us suppose that (xx,X2) <E ^(2). Some technics already used in 
the previous section give: 



(31) 



At this point of the proof, two cases must be distinguished: either both 
local configurations r\\ and 7/2 are equal to 77 or not. In the first case, 
Lemma \l .41 and (131(1 imply: 



-2d\2 
I A 



Then, the quantity E 2 (rj, 77) which is actually the second moment of X n {rj), 
is bounded by A 2 . In the other case, at least one of the two local configu- 
rations jji, r\% 6 D r {r\) is different from 77, i.e. has at least k(j)) + 1 positive 
vertices. Then, coupling Lemma \'2 . 21 with (|31|l . it follows that: 



E 2 (m,V2) < Mle 2a{n){k{71l)+k{7l2) ~ 2k{7l)) 
which is a 0(e 2Q ( n >) = 0{nr d / ki -^). 



□ 
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